Factorization fits to charmless strangeless B decays 
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We present fits to charmless strangeless hadronic B decay data for mean branching ratios and 
CP-violating asymmetries using the QCD factorization model of Beneke et al. Apart from one 
CP-violating parameter, the model gives a very good representation of 26 measured data. We find 
the CKM angle a = (93.5 ± 8.4 — 1.3)° and to be quite stable to plausible "charming penguin" 
corrections. 
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I. INTRODUCTION 

A wealth of experimental data on hadronic charm- 
less B decays has become available from the BaBar and 
Belle experiments. Many new branching ratios and CP- 
violating parameters have been measured within more 
precise error limits. These studies of the numerous B de- 
cay channels are designed to test the Cabbibo-Kobayashi- 
Maskawa (CKM) explanation of CP violation in the stan- 
dard model. 

We have made previous attempts 0, to understand 
charmless B meson decay data in terms of the QCD fac- 
torization model of Beneke et al. (BBNS) U 3- In 
fitting the data we found evidence for charming penguin- 
like contributions Q to the decay amplitudes in addition 
to the BBNS amplitudes. The charming penguin contri- 
bution to the strangeless modes is suppressed relative to 
the strange modes by one power of the Wolfenstein pa- 
rameter A K, 0.22. In P] we attempted a simultaneous fit 
to both strange and non-strange channels and, although 
we obtained a satisfactory fit to measured branching ra- 
tios and some CP- violating asymmetries, the predicted 
CKM angles deviated significantly from other analyses. 
The strange channels in isolation provide the best data 
for examining the phenomenology of charming penguins 
1^ . With this information their influence on the strange- 
less modes should then be investigated. 

We report here the results of our application of a BBNS 
analysis to charmless strangeless B decays. We use the 
method and notation of and attempt to fit the data 
on the mean branching ratios and CP- violating asymme- 
tries of the decays B — > tttt, 7rp, ttw and pp, together with 
sin(2/3), 26 independent measurements in all. We include 
many more CP asymmetry measurements in our data 
set and also three new pp channels which exhibit the ex- 
pected helicity-zero dominance. The paper is organized 
as follows. In Sec. II we briefly discuss the structure 
of the decay amplitude within QCD factorization. The 
weak annihilation contributions are summarized in Sec. 
Ill and the method and results of our best fit are pre- 
sented in Sec. IV. Sec. V contains our discussion and 
conclusions. 



II. DECAY AMPLITUDE IN QCD 
FACTORIZATION 

In QCD factorization, the matrix elements of the op- 
erators in the effective Hamiltonian Ticff are separated 
into short distance contributions at scale 0(l/TOf,) that 
are perturbatively calculable and long distance contribu- 
tions 0(1/Aqcd) that are parametrized. The amplitude 
for B decay into two light hadrons (mesons) Afi2 has the 
form, neglecting weak annihilation processes, [EH 111 

{M,Kh\n.s\B) = ^{[E ^««"+ E E 

^ i=l,2 p=ti,ci=3,...,6,9 

X [r,(Afi, M2) + T,(M2, Ml)]} (1) 

where Ap = VpdYpb is a product of CKM matrix elements. 
In Eq. ^ we have included current-current tree pro- 
cesses represented by the Wilson coefficients Ci,2, QCD 
penguin processes represented by 0^^,,,^ and the domi- 
nant electroweak penguin process represented by Cg . The 
factorization matrix elements Ti{Mi, M2) involve prod- 
ucts of two-quark current matrix elements and, neglect- 
ing factors {mp,v /nT-B)'^ and assuming zero helicity dom- 
inance for B ^ V1V2 decays, have the form ^] 

T,(Mi, M2) = c m%fM,FM, (2) 

where /m are the well determined electroweak decay con- 
stants {fjr, fuj, fp}, Pa/ are the B transition form factors 
{F,^, A^, Ap} and the constant c is a product of factors 
1, ±1/V2, e.t.c. from the flavor composition of the B and 
Ml, 2 mesons. 

The factorization coefficients Oi are calculated from the 
Wilson coefficients Ci{p,) at a scale p of 0(mb) and have 
the form |3] 

a,(Mi,M2) = a^,/(M2) -l-a,,/7(Mi,M2) (3) 

where Mi is the recoil meson containing the spectator 
(anti) quark and M2 is the emitted meson. The complex 
quantities aij describe the formation of M2, including 
nonfactorizable corrections from hard gluon exchange or 
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light quark loops in penguins. The hard gluon exchanges 
with the spectator quark are described by the (possibly) 
complex quantities aijj. In these correction terms the 
leading-twist light cone distribution functions for both 
pseudoscalar and vector mesons are expanded in the first 
few terms of a Gegenbauer expansion 

<fMix, ^l) = 6x(l - x)[l + Y.^'JCl'\2x - 1)]. (4) 

n 

The light cone corrections are anticipated to be small 
but they are not well established. Consequently it is 
common to use the asymptotic form <^m{x, ^) = 6x(l — 
x), valid for the mass scale /x — > oo, in applications of 
QCD factorization. All coefficients aij except ag,/ are 
then the same for all decays. If the corrections af^^ ^'^^ 
included, the coefficients most affected are a2j and a^j 
(see Table IV of Q). 

The contributions of the a^j/ coefficients to the decay 
amplitudes have the form 



47r 
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where i' = i — (—1)*, £i = 
— 1, eg = 0, and 
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Here r*^^ are chiral factors, Xf^ = Xh and 0| 



III. ANNIHILATION CONTRIBUTIONS 

B meson decay can also be initiated by b quark an- 
nihilation with its partner. Although the annihilation 
contributions to the decay amplitude are formally of 
0(AQCD/wb) and power suppressed, they violate QCD 
factorization because of end point divergences. However 
these weak annihilation contributions can be included 
into the decay amplitudes by treating the end point di- 
vergences as phenomenological complex parameters Xa- 

The annihilation contribution to the decay amplitude 

is 



{NhM,\-HlT\B) 



where 



Cf 

BmiM2 = J^fBfMjM2, 



(10) 



(11) 



Cp — {N^ — l)/2Nc and Nc is the number of colors. 
The quantities ^^^'{(Mi, M2), where the superscript i{f) 
denotes gluon emission from initial (final) state quarks, 
evaluated using the asymptotic form of Eq. Q are given 
by i 



"iiaX + aX)XH - (6 + 9af + lla^ 



(7) 



The non-perturbative complex parameter Xh is the con- 
tribution of a logarithmic end-point divergence in the in- 
tegration over the twist-3 light cone distribution function 



Xh = 



dx 



1 



and is usually parameterized as 



XH = \n 



rriB 



A 



QCD 



iipH 



(8) 



(9) 



where ln(mB/AQCD) — 3.03. The coefficients Uiji are 
not universal even when light cone corrections are ne- 
glected and they contain the parameter Xh which is 
only loosely constrained by model estimations. These 
Qiji contributions to the decay amplitudes are indepen- 
dent of the B transition form factors but do involve the 
poorly determined parameter Jb/^b where Jb is the B 
leptonic decay constant and \b ~ 0.35 GeV is related 
to the B light cone distribution function. We note that 
substantial light cone corrections af^2 significantly 
enhance the aiji coefficients. 

The energies involved in the calculation of aiji im- 
ply that the appropriate scale is not that of the scale 
u used in calculating the a^j but — V^ftM where 
[3 A/j = 0.5 GeV. For the choice /i = m}j/2 this gives 
asfBl{mB\B) = 0.0209. 
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{X\ - 2Xa)\ , 



^ ){X\-2Xa^-\ 



4(^^1.1^2) - 18™,(r];i +r^^)(2XA-l)(XA-2)(12) 



The coefficients di{M\, M2) are Clebsch-Gordan type fac- 
tors and are tabulated in [J for all the decays studied 
here apart from B P~ for which the required d- 

coefficients are (1, 0, —1, 2, —1, 2). 
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IV. FITTING METHOD 

We attempt to fit the experimental data on the mean 
branching ratios and CP- violating asymmetries of the 
decays B tttt, 7rp, ttw and pp, together with sin(2/3). 
This gives 26 independent measurements in all. There are 
many parameters in the equations of the BBNS theory 
that are not precisely known, we choose the B semilep- 
tonic transition form factors F^r , Ap and A,^ , the Wolfen- 
stein CKM parameters p and 77 and the complex non- 
perturbative hard scattering parameter and annihi- 
lation parameter Xj^ as our nine fitting parameters. 

We assign to each of the 26 independent data a number 
a and construct a function of the nine fitting param- 
eters 

26 

= ^[((Theory)^ — (mean data value)Q) /ctq]^. (13) 

ctq, is an experimental error formed by amalgamating sta- 
tistical and systematic errors. We then use a minimiza- 
tion procedure based upon the program MRQMIN from 
to try to obtain acceptable fits of the theory to ex- 
periment. An acceptable fit has a low value of a-t the 
minimum with parameter values that lie within known 
acceptable limits. Other parameters in the analysis such 
as the correction factors to the light cone distribution 
functions of the participating mesons are poorly known 
and not easy to incorporate as variables in the minimiza- 
tion procedure. In this paper we use three sets of values 
for af^2i the zero set and two sets suggested by theory 
(those of Q and Q) to examine the sensitivity of the 
results to the choice of these correction factors. We also 
examine how the results depend on our choice of the chi- 
ral factor and the factor = otsfs / {"niB^B) which 
enters into the coefficients aiji representing the specta- 
tor quark interactions. 

Table shows the best fit values for the parameters, 
together with the statistical error from the and a sys- 
tematic error from the different light cone correction fac- 
tors, e.t.c. This fit has a minimum of 9.2 and is for the 
light cone distribution functions of Lii and Yang 8] , the 
chiral factors = 1.0 and = 0.1, and = 0.0314. 
An unexpected feature of this fit is the large value of 
the spectator quark interaction term pu which was not 
expected to be greater than 3.0. It can be seen how- 
ever that its statistical error is large, the function 
has a very shallow worm hole through parameter space. 
Constraining pu to be 3.0 gives a minimum x^ of 12.5 
and a quite acceptable fit. The parameter changes this 
induces, the changes due to the different light cone dis- 
tribution functions, and those from taking r'^ = 0.9 or 
(^B = 0.0209, are summarized in the column of system- 
atic errors. These errors, like the statistical errors, are 
of course highly correlated but they indicate the sensi- 
tivity of the parameters to different choices. Our fit also 
implies that Vub — (3.89 ± 0.24) x lO"-^, to be compared 
with (3.7 ±0.8) X IQ'^. 



TABLE I: Best fit values and errors for the B transition 
form factors Ft,, Ap and A^, the complex hard scattering 
Xh ~ 3.03 -I- ph exp{i(j)H) and annihilation Xa — xa + iyA 
parameters, and the CKM angles acKM and /3ckm. 



Parameter 


Mean 


Statistical 


Systematic 




VclIug 


Error 


Error 




0.218 


±0.022 


+0.015, -0.003 


Ap 


0.317 


±0.027 


-0.021 




0.372 


±0.050 


-0.03 


PH 


8.73 


±2.82 


-^3.1,-5.26 


4<H 


-2.05 


±0.15 


-H0.49, -0.16 


XA 


-1.01 


±0.50 


-h0.02, -0.34 


VA 


2.89 


±0.85 


-hO.73, -0.39 


CfCKM 


93.5° 


±8.4° 


-1.3° 


/3CKM 


24.2° 


±2.3° 


+0.2°, -0.4° 



The predicted mean branching ratios and CP asym- 
metries for the parameter values of Table are shown 
in Table (nj together with the individual contributions 
to the total minimum x^ of 9.2. Also shown are the es- 
timated theoretical errors, the statistical errors are the 
Icr standard deviations from the 9x9 error matrix of 
the x^ minimization and the systematic errors arise from 
changes in the light cone distribution functions, e.t.c. 

If significant contributions, such as charming penguins, 
are needed to understand the charmless strange B decays 
then it is to be expected that similar but smaller contri- 
butions should be added to the BBNS analysis. In the 
spirit of (li we have investigated the influence of charming 
penguins, as will be discussed in the next section. 

V. DISCUSSION AND CONCLUSIONS 

Apart from the value of pu , which has been discussed 
above, the fitting parameters of Tableware within their 
expected ranges. In particular, the CKM parameters V^b 
and the angles acKM and /3ckm have very acceptable 
values within quite tight bounds. All the experimental 
data, apart from C^rTpi are well fitted by the analysis. 
The largest individual x^ of 3.76 is for this CP-violating 
parameter and it can be seen that our analysis gives 
a value much smaller than the measured value, which 
is from the BaBar experiment. The theoretical result 
is in much better accord with the Belle measurement 
C^Tpi = 0.25 ± 0.17+° °g 10] (talk by M. Giorgi). A 
remeasurement of this parameter would be of interest. 

The angles acKM and /3ckm are the most important 
parameters to come out of the analysis. It is of in- 
terest to compare our value of acKM with those from 
isospin analyses and the Grossman- Quinn bound. These 
analyses are independent of the details of the QCD 
penguin contributions. The small p'^ p'^ branching frac- 
tion makes the pp analysis the most precise and yields 
acKM = 96° ± 10° ± 4° ± 11° [13 (talk by M. Giorgi) 
where the last error is from the Grossman- Quinn bound. 
However our neglect in particular of charming penguin 
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TABLE II: Measured branching ratios and CP asymmetries (Data), experimental error (a), best fit theoretical values (Theory), 
estimated theoretical errors (Statistical) and (Systematic), contribution Xa to the total minimum of 9.2, and the reference 
(Ref.) for the experimental data. All branching ratios are in units of 10~^. 





Data 


O" 


Theory 


Statistical 
Error 


Systematic 
Error 




Ref. 


Br(7r+7r") 


4.6 


0.4 


4.73 


±0.45 


+0.22,-0.07 


0.11 


8 


Br(7r°7r°) 


1.2 


0.4 


1.15 


±0.33 


±0.0, -0.38 


0.02 


9 


Br{uj7v°) 


< 1.0 




0.48 


±0.30 


±0.02,-0.23 


0.27 


8 


Br(7r±p=F) 


24.0 


2.5 


23.6 


±1.54 


±0.20,-0.06 


0.03 


8 


Br(7r°p°) 


5.1 


2.4 


2.90 


±0.44 


±0.0,-1.08 


0.84 


8 


Br(p+p-) 


30.0 


6.0 


31.9 


±4.02 


±0.0, -1.90 


0.10 


9 


Br(p'V") 


< 1 




0.66 


±0.23 


±0.0, -0.33 


0.12 


9 


Br(7r±7r°) 


5.2 


0.8 


4.51 


±0.65 


±0.71,-0.24 


0.74 


8 


Br(tJ7r±) 


5.9 


0.8 


5.90 


±1.00 


±0.0 


0.00 


8 


Br(7r±p'') 


9.2 


1.2 


10.1 


±1.03 


±0.40, -0.0 


0.61 


10 


Br(7r°p±) 


12.0 


1.9 


10.8 


±0.80 


±0.30,-1.70 


0.42 


8 


Br(p°p±) 


26.4 


6.4 


23.0 


±2.19 


±0.0,-3.10 


0.28 


11 


Ctt+tt- 


-0.09 


0.19 


-0.13 


±0.03 


±0.01 


0.04 


9 


s + 


-0.30 


0.21 


-0.21 


±0.14 


±0.0,-0.06 


0.17 


9 




0.05 


0.15 


-0.02 


±0.0 


±0.01,-0.0 


0.23 


9 




-0.17 


0.13 


-0.15 


±0.08 


±0.0,-0.03 


0.03 


10 


AnOp± 


0.23 


0.22 


0.15 


±0.06 


±0.03, -0.0 


0.13 


10 






U. ID 


u.uo 


-1-0 m 




0.(0 


Q 




0.15 


0.14 


0.18 


±0.08 


±0.0,-0.01 


0.25 


9 


'S',rTp± 


-0.10 


0.18 


-0.15 


±0.15 


±0.01 


0.08 


9 


AS'^Tp± 


0.22 


0.18 


0.15 


±0.05 


+0.02,-0.11 


0.15 


9 




-0.09 


0.06 


-0.10 


±0.03 


±0.0 


0.04 


9 


CpTpi 


-0.23 


0.38 


0.03 


±0.09 


±0.0 


0.47 


9 


'S'pTpi 


-0.19 


0.44 


-0.27 


±0.02 


+0.04, -0.0 


0.03 


9 


^pTpO 

sin(2/3) 


-0.09 


0.16 


-0.01 


±0.0 


±0.0 


0.27 


11 


0.734 


0.054 


0.747 


±0.046 


±0.005 


0.06 
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contributions to the decay amplitudes will bias our es- 
timates. We have only made a preliminary assessment 
of the influence of charming penguins, our conclusions so 
far are that the small branching fractions like Br(7r°7r°) 
are somewhat sensitive but the CP-violating data can 
be more affected, especially those like C^Tp± which our 
analysis predicts to be very small. So far we have not 
found the angles acKM and /3ckm to be particularly sen- 



sitive to charming penguins but more work is required to 
get reliable estimates of the systematic errors. 

In conclusion, we find that the BBNS analysis so far 
does well in fitting the data. We see no evidence in 
charmless strangeless decays for any physics beyond the 
Standard Model and the analysis holds the promise of 
giving tight constraints on the values of the CKM pa- 
rameters. 
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